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Abstract 

A monochromatic line in the cosmic neutrino spectrum would be a smoking 
gun signature of dark matter. It is intriguing that the IceCube experiment has 
recently reported two PeV neutrino events with energies that may be equal 
up to experimental uncertainties, and which have a probability of being a 
background fluctuation estimated to be less than a percent. Here we explore 
prospects for these events to be the first indication of a monochromatic line 
signal from dark matter. While measurable annihilation signatures would 
seem to be impossible at such energies, we discuss the dark matter quantum 
numbers, effective operators, and lifetimes which could lead to an appropriate 
signal from dark matter decays. We will show that the set of possible decay 
operators is rather constrained, and will focus on several viable candidates 
which could explain the IceCube events; R-parity violating gravitinos, hidden 
sector gauge bosons, and singlet fermions in an extra dimension. In essentially 
all cases we find that a PeV neutrino line signal from dark matter would be 
accompanied by a potentially observable continuum spectrum of neutrinos 
rising towards lower energies. 



1 Introduction 



The IceCube collaboration has very recently reported a detection of two neutrino 
events with energies of 1.1 PeV and 1.3 PeV in an energy range where no more 
than 0.01 background events was expected from atmospheric neutrinos |U [21 E]. 
These are stated to be either electron neutrino charged current events, or neutral 
current events of any neutrino flavor. It is interesting that the two detected neutrinos 
have such similar energies, and indeed, most astrophysical sources are expected to 
produce power-law spectra; one might expect to see more events at lower energies, 
above the experimental threshold of 0.4 PeV, as well as at higher energies, around 
6.8 PeV, where the detector sensitivity is enhanced by the Glashow resonance [I]. 
The data may thus suggest a peak in the neutrino spectrum around 1 PeV. It is 
possible that such a spectrum could be produced by some astrophysical sources [5] 
[TJ, including intergalactic interactions of cosmic rays produced by blazars [5] [16], 
but these models rely on some assumptions about the properties and evolution of 
the sources, as well as the intergalactic magnetic fields. 

The IceCube observations raise a question of whether dark matter could be com- 
posed of relic particles whose decays or annihilations into neutrinos produce a line- 
like feature in the neutrino spectrum at ~1.2 PeV energy. Similar to a line in the 
gamma ray spectrum, a line in the neutrino spectrum could be considered a "smoking 
gun" signature for dark matter. Such line-like neutrino signatures from dark matter 
have been considered before [T7J-[TS], but in this paper we concentrate on the pos- 
sibility of obtaining such a signal at the PeV scale, where the dark matter particle 
cannot be a simple thermal relic. As we will show, the possibilities for obtaining a 
neutrino line signal from dark matter at such energies are highly constrained, but 
there are nevertheless various viable scenarios. 

For dark matter with an annihilation cross section into neutrinos saturating the 
unitarity limit, o"Ann < 47r/(mp M w 2 ), the event rate expected at a neutrino telescope 
of fiducial volume V and nucleon number density tt-n is of order 

TEvcnts ~ V L MW n N cr N ( ] (a Ann y) < 1 per few hundred years, (1) 

\ m DM / 

where we have taken the neutrino-nucleon scattering cross section to be <7n ~ 9 x 
10 _34 cm 2 at E u ~ 1.2 PeV [20], and the nucleon number density to be that of ice, 
tin — ^icc — 5 x 10 23 /cm 3 . Pdm-i v, and Lmw are the milky way dark matter density 
(taken near the Earth for the purpose of our estimate), the typical dark matter 
particle velocity, and the rough linear dimension of our galaxy, where these are fixed 
to be 0.4 GeV/cm 3 , 10 -3 , and 10 kpc, respectively. The fiducial volume V is set to 
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be 1 km 3 , which is roughly the size of the IceCube detector. We see that obtaining 
a neutrino line signal from dark matter annihilations at the PeV scale is essentially 
not possible. In what follows we will therefore restrict ourselves to the possibility of 
a signal from dark matter decays. 

For dark matter decays also, obtaining a neutrino line signal at the energies of 
interest here turns out to be challenging. Indeed, suppose one wishes to mediate an 
appropriate decay via a simple dimension 4 operator such as C D XipLH, where A is 
a coupling constant, ip is the dark matter particle, L is a lepton doublet, and H is 
the Higgs doublet. Then the decay rate to neutrinos is Tdm = j^t^bm- Similarly to 
the annihilation case above, we may estimate the event rate at a neutrino detector 
for rriDM — 1.2 PeV as 



Clearly an exceptionally tiny coupling is required to obtain an appropriate signal, 
and a certain amount of model building would appear necessary. 

We may also consider whether or not higher dimension operators, suppressed 
by some large mass scale, could give more naturally small event rates. For higher 
dimension operators, however, it is a nontrivial constraint that in order to obtain 
a line signal, the decay final state must be two-body. Indeed, for many interesting 
operators, neutrinos appear in the gauge singlet combination LH, and although 
naively this could lead to a neutrino decay with H replaced by its vacuum expectation 
value v, this tends not to be the dominant process due to the large dark matter masses 
under consideration. For example, if one considers the operator £ D <ft(LH) 2 /A for 
a scalar dark matter particle 0, and with A a heavy mass scale, then the square of 
the amplitude for a four-body decay with two neutrinos and two Higgses is larger 
than that for a two-body neutrino decay by a factor of ~ (^dmAO 4 - For heavy dark 
matter masses, phase space suppressions for multi-body final states are not enough 
to prevent the four-body decay from being by far dominant. 

In this paper we will comprehensively discuss effective operators which could 
mediate the decays of heavy dark matter particles into monochromatic neutrino lines, 
and we will find that only a handful of operators are viable. Several of these stand 
out as being particularly interesting, and we will discuss possible models for them in 
detail. These will include the cases of gravitino dark matter, with a mass motivated 
by the recent 125 GeV Higgs discovery, a hidden gauge boson with an extremely 
small mixing with hypercharge, and a singlet fermion in an extra dimension. In 
each case we will discuss simple ways in which an appropriately long lifetime for 



P Events ~ V LmW iN^N 



Pdm 




(2) 
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the dark matter particle may be obtained in order to explain the IceCube data. In 
the gravitino case, the decay operator may be naturally suppressed by the scales 
of R-parity violation and lepton number violation. In the gauge boson case, the 
kinetic mixing with hypercharge may be suppressed by the scale of non-abelian 
gauge symmetry breaking in the hidden sector, as well as the breaking of grand 
unified symmetry in the visible sector. In the extra dimensional model, the required 
highly suppressed coupling may be produced naturally by an exponentially small 
wave-function factor. 

The outline of this paper is as follows: In section [2] we will review the nature of 
the PeV IceCube neutrino events, as well as discuss the lifetime and mass of dark 
matter particles which may be able to explain them. In section [3] we will discuss in 
general the effective operators which might be able to lead to an appropriate dark 
matter decay. Models yielding some of these operators will be discussed in section 
IH An interesting conclusion of our analysis will be that in essentially all cases, a 
monochromatic neutrino line would be accompanied by a continuum spectrum of 
neutrinos rising towards lower energies. The prospects for detecting such a signature 
will be discussed in section |5j 

2 IceCube Events 

Before going on to discuss effective operators which could mediate the decays of 
heavy dark matter particles into monochromatic neutrino lines, we summarize the 
situation with the PeV neutrino events which have recently been reported by the 
IceCube collaboration. According to a plot in reference j2], the exposures at the 
energies of the two events turn out to be 4.4 x 10 9 [m 2 s sr] and 5.9 x 10 9 [m 2 s sr] 
for the 1.1 and 1.3 PeV events, respectively, assuming that both the two events were 
caused by electron neutrinosr 1 ! It follows that the total flux may be estimated to 
be 

T ~ 4.0 x 1(T 14 [cm" 2 s" 1 sr^ 1 ]. (3) 

The observed neutrino event energies imply a mass for the dark matter particle of 
about 2.4 PeV, while the neutrino flux can be related to the lifetime for dark matter 
neutrino decays. When the mass of the decaying dark matter particle is assumed to 

^ tl It is possible that one or both events could have been caused by neutral current interactions 
of arbitrary flavor, but in such cases one would expect the event energies to be much more spread 
out. 
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be 2.4 PeV, the predicted flux of line neutrinos is estimated to be 

dT /in 29 s\ 

El-!- ~ 9.5 x lO" 3 ^ 8{E V - m DM /2) [GeV cm" 2 s" 1 sr^ 1 ], (4) 

dhj v \ r DM / 

where E u , imdm, tdm ; and N u are the neutrino energy, the dark matter mass, its 

lifetime, and the number of neutrinos produced in each decay, respectively. Here 

the NFW profile was used for the dark matter density in our galaxy, and we have 

adopted profile parameters with a critical radius of r c = 20 kpc, and a density at 

the solar-system of p Q = 0.39 GeV/cm 3 [21]. The total flux is then given by 

/10 29 s\ 

T ~ 0.76^ x 10~ 14 — — [cm" 2 s" 1 sr- 1 ]. (5) 

V r DM / 

By comparing this prediction with the flux in equation flS]), we find that the the 
lifetime of the dark matter particle must have the following value in order to explain 
the data: 

r DM ^ 1.9JV„ x 10 28 s. (6) 

We have thus found that a decaying dark matter particle with a mass of about 
2.4 PeV and with a lifetime as given in equation ([6]) can explain the IceCube PeV 
neutrino events. Note that as a result of neutrino oscillations, all neutrino flavors 
will contribute equally to the final signal, independent of the original flavor structure 
of the dark matter decays. 



3 Effective Operators 

Here we list all operators which might lead to a high energy monochromatic neutrino 
line from dark matter decays. In table [U we show possible dark matter candidates, 
defined by standard model SU(2)i and U(l)y quantum numbers. We only list candi- 
dates which have a leading decay operator to two standard model particles, including 
at least one neutrino. We exclude cases in which there is an alternate decay mode 



through an operator of lower dimensionP 2 ! or which require decays to additional 



non-standard model particles. Note that case 1 is a slight exception to this rule, 
since a decay through a lower dimension operator H^<pH c is possible, but we include 
this case since the two Higgs decay may be forbidden by lepton number. In the final 
column of the table, we give the coefficient for the operator required to explain the 
two IceCube events based on the flux in the previous section. 



# 2 In cases in which the dark matter particle carries hypercharge (case 1, 4, and 5 in table [1]), a 
Dirac mass partner is required. We only include in the table operators of lowest dimension when 
considering all operators allowed for cither member of the Dirac pair. 
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Case 


Spin 


SU(2) L 


U(l)y 


Decay Operator 


Coefficient for IceCube Data 


1. 





3 


1 


L c (pL 


9.5 x 10" 30 


2. 


1/2 








LH c ip 


2.7 x 1(T 29 


3. 


1/2 


3 







3.8 x 10~ 29 


4. 


1/2 


2 


-1/2 


LFip 


5.6 x 1(T 30 (PeV" 1 ) 


5. 


1/2 


3 


-1 


LiJj a T a H 


2.7 x 1(T 29 


6. 


1 








LfL 


3.3 x 1(T 29 


7. 


3/2 










1.9 x 10~ 29 (PeV- 1 ) 



Table 1: Daric matter candidates and the decay operators that may lead to a mono- 
chromatic neutrino line signature. Here, L and H represent the SM lepton and Higgs 
doublets, respectively, while the dark matter particle is labeled by (ft, ip, or tp^ 1 , de- 
pending on whether it has spin 0, 1/2, 1, or 3/2. The notation F in case 4 denotes either 
B^o^ , B^cr^, W^ u T a cr^ or W£ u T a cr^ w ith B^ and W* v being the held strength ten- 
sors of the SM U(l )y and SU(2)l gauge Eelds. In the hnal column, we give the coefficient 
for the operator required in order to explain the two anomalous neutrino events reported 
by the IceCube collaboration, assuming a dark matter particle mass of 2 A x 10 6 GeV. 

Cases in which the dark matter particle carries electric or color charge have not 
been included in the table. Electrically charged dark matter is severely constrained 
by several observations and experiments, and is required to be heavier than about 
10 12 GeV [221 1231 121] , primarily by difficulties with structure formation. Colored dark 
matter, similarly, must be heaver than about 10 16 GeV [251 126] . with the primary 
constraint coming from the possibility of overheating the Earth's core. We have, 
on the other hand, included cases with non-zero hypercharge, which naively have 
excluded tree level Z-boson exchange signatures at dark matter direct detection ex- 
periments. These constraints can be avoided, however, if there is a higher-dimension 
operator which induces a splitting among the components of the dark matter field 
in such a way that the lightest state becomes a Majorana particle. Such a splitting 
is then required to be larger than the recoil energies produced at direct detection 
experiments. This is in fact what occurs for the case of higgsino dark matter in 
supersymmetric models- the mixing with Majorana gauginos causes the splitting. 

As discussed in the introduction, all decay operators we consider in table[T]contain 
only three fields. This was done in order to ensure that a monochromatic neutrino 
line signal dominates over other decay modes. Operators requiring extra insertions 
of Higgs vacuum expectation values to yield a monochromatic neutrino decay are 
not allowed, since multi-body decays with extra Higgs particles would give overly 
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large alternate cosmic ray signatures. 



4 Models 

4.1 Gravitino Dark Matter with R-Parity Violation 

Our first example model comes from the operator listed clS CclSG 7 in table [TJ This 
operator requires the dark matter particle to have spin 3/2- namely, to be a gravitino. 
Here we will show that, in an R-parity violating context, it is straightforward to 
obtain a monochromatic neutrino line from gravitino dark matter decays, with a 
mass and lifetime appropriate for explaining the PeV IceCube events. 

We begin by considering the mass of the gravitino. Both the ATLAS and CMS 
collaborations at the LHC have reported the discovery of a Higgs boson with mass 
around 126 GeV [2SII2S]. The mass is somewhat heavier than one could expect in the 
minimal supersymmetric standard model (MSSM), but large radiative corrections to 
the Higgs quartic coupling [30]- [33] can lead to a heavier lightest Higgs mass in the 



MSSMF 3 !. When the left-right mixing of scalar top quarks is negligible and tan /3 — 2, 
the typical scale of sparticle masses must be 0(1) PeV assuming simple gravity 
mediated SUSY breaking. Based on this observation, several concrete models have 
been proposed [35]- [37], which are attractive from the viewpoint of the SUSY-Flavor 
and CP problems because all dangerous flavor changing processes are suppressed by 
heavy sfermion masses. Gravitino dark matter with a PeV mass is, therefore, quite 
consistent with the observed Higgs mass under the assumption that the gravitino is 
the lightest supersymmetric particle (LSP). 

Let us now consider the lifetime of the gravitino, whose decay must be induced by 
some form of R-parity violation. Here we will consider a simple set of assumptions 
which will imply that the leading R-parity violating operator in the superpotential 
will be of the form LH U , with a coefficient of order m^ 2 /M p \, where L is a lepton 
doublet of arbitrary flavor, m 3 / 2 is the gravitino mass, and M p \ is the Planck scale. 
There are several assumptions required: The first is that the R-charges of all MSSM 
matter fields are equal to 1, while those of the Higgses are equal to 0. Next we sup- 
pose that it is actually a Z3 subgroup of U(1)r which is a symmetry of the theory, and 
not the full continuous U(1)r. Since the gravitino mass is a spurion for R-symmetry 
breaking with R-charge 2, and Z 3 R-symmetry requires that superpotential terms 
have R-charge equal to 2 mod 3, we find that the R-parity violating operator LH U 



* 3 It is also possible that large supersymmetry breaking terms cause some squarks to form Higgs- 
like bound states, hence relaxing the MSSM limits on the mass of the lightest Higgs boson [33]. 
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appears with a coefficient of m\, 2 /M p \ as promised [38]. Note that other R-parity vi- 
olating operators, UDD, LLE and QLD all appear at order m 3 / 2 /M p i, and therefore 
also with one suppression by the Planck scale. These will lead to continuum neutrino 
decay spectra in addition to the monochromatic line (plus continuum) obtained from 
the R-parity violating operator LH U . Note, however that the decay rates from these 
other R-parity violating operators will have additional phase space suppressions due 
to extra final state particles, and are thus naively expected to be sub-dominant. 

As a result of the LH U operator, the lifetime of the gravitino to decay into a 
neutrino plus a Higgs or a neutrino plus a Z boson is estimated to be [39] 

r 3/2 ~ 192tt (M pl /m 3/2 ) 4 ~ 10 20 s, (7) 

where for illustration the sneutrino mass is assumed to be the same as the gravitino 
mass, namely 2.4 PeV, and M pl ~ 2.4 x 10 18 GeV. This lifetime of order 10 20 s is 
too short to be consistent with the two IceCube PeV neutrino events. 

We now note, however, the discussion leading to equation (J7|) potentially misses 
an important point. Indeed, the operator LH U carries B— L charge —1. In any 
case, B— L symmetry must be broken to allow for Majorana masses for right handed 
neutrinos Nr [10]- [42]. Therefore, as is standard we may introduce B and B fields 
carrying B— L charges +1 and —1, respectively in order to break this symmetry. Now 
one can write W D NrNr(B) 2 / M p \, which provides the Majorana mass term. If one 
assumes M N ~ 10 10 GeV, the expectation value of the B field is (B) ~ 10~ 4 M pl . 



Then the R-parity violating operator becomes of ordeir 4 

W={ml /2 (B)/Ml)LH u . (8) 

Because of this modification, the lifetime of the gravitino (decaying into a neutrino 
and either a Higgs boson or Z boson) is now about 10 8 times longer than the lifetime 
in equation (17j) . To be more precise, the lifetime is then estimated to be 

r 3/2 ~ 1927r(A.f pl /( J B}) 2 (M pl /m 3/2 ) 4 m- / 1 2 ~ 10 28 s, (9) 

where the Majorana mass has been set to 10 10 GeV. This lifetime is fully consistent 
with the one implied by the IceCube PeV neutrino flux in equation ([6]). 

Finally, let us discuss gravitino production in the early universe and the dark 
matter abundance. The next-to-lightest-superpartner (NLSP) in this model will 
generically decay to the gravitino (plus its standard model partner) with a lifetime 
of order M^/m^^p. With an NLSP mass at the PeV scale this is roughly of order 



# 4 Note that B and B have R charge according to our charge assignments, as required. 
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10 -5 seconds. The NLSP decays thus do not disrupt the successful predictions of big 
bang nucleosynthesis. On the other hand, the NLSP freeze-out abundance, which 
will then be converted into the gravitino relic abundance, will be too high by perhaps 
a factor of ~ 10 5 due to the large NLSP mass. If the reheating temperature is above 
the NLSP mass, entropy production by a factor of ~ 10 5 will thus be required to 
dilute the dark matter abundance. If the reheating temperature is below the NLSP 
mass on the other hand, then it is possible to produce an appropriate gravitino 
abundance through a small branching fraction of the inflaton into the gravitino. 



4.2 Hidden Sector Gauge Boson 

Another interesting possibility for a dark matter particle which could give a neutrino 
line signature at IceCube comes from case 6 in table [TJ Here we require a new gauge 
boson V^ 1 with a very small coupling ~ 10~ 28 to at least one standard model lepton. 
What is very interesting about this case is that a coupling of this size may be 
obtained in a very simple and natural wayP 5 ! In particular, let us suppose that 
the visible sector is part of a standard grand unified theory, with a unification scale 
of Mqut ~ 10 16 GeV. We may take a minimal SU(5) theory with GUT symmetry 
broken by the vacuum expectation value of an adjoint scalar field S for illustration. 

Now, we consider the possibility that there is a completely hidden sector with a 
new non-abelian gauge symmetry, broken at the PeV scale. Let us take this gauge 
symmetry to be SU(2) for simplicity, and suppose that it is completely broken by the 
vacuum expectation value (vev) of a scalar field <3> in the fundamental representation. 
Because both the visible and hidden sector gauge symmetries are fundamentally 
non-abelian, dimension four kinetic mixing between their respective field strengths 
coming from an operator ~ F fJLU V^ u is forbidden, where now and V^ u are taken 
to be the SU(5) and hidden sector field strengths, respectively. However, after gauge 
symmetry breaking, such mixing is induced by Planck suppressed operators, even 
if there are no new particles carrying both visible and hidden quantum numbers. 
In the present example, the minimal Planck suppressed operator which results in 



mixing between the visible and hidden gauge bosons is given 



£ D Jp^F^&V^®. (10) 

This leads to a kinetic mixing between hypercharge and the lightest new gauge boson 

of order ^J s ^ ~ PeV ,^f GUT ~ 10~ 28 . The hidden gauge boson will then obtain a 

P i P i 



# 5 For another model which may be used to give a similar resulting decay operator please see 
^Operators with different combinations of $ and $t arc similarly allowed. 
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coupling to the standard model leptons of the right order to explain the IceCube 
data! Of course, we are assuming here that there are no light hidden sector particles 
into which the hidden gauge boson may rapidly decay. 

Note that in the absence of supersymmetry, this model introduces a new hierarchy 
problem for the mass of the scalar field $. There is also a dangerous allowed quartic 
coupling between $ and the visible sector Higgs boson, which would lead to very 
rapid decays to Higgs bosons. There is, however, no obstacle to implementing 
the model in a supersymmetric framework, and doing so can prevent the quadratic 
divergence of the $ mass, as well as forbid the $/Higgs quartic coupling. On the 
other hand, there is one additional type of dangerous operator which supersymmetry 
cannot forbid. This is an operator of the form 

^V^H%H, (11) 

which may be generated by a Kahler potential term -rj T &$H^H and which results in 

Ju pi 

decays to two Higgses. Note that the operator (fTUj) which leads to the monochro- 
matic neutrino line is suppressed by an additional factor of J g^ T compared to (fTTT) . 
We thus require that the new operator be suppressed by a factor of about 100-1000 
beyond the naive estimate in ( fTTT) of jjt in order that the monochromatic neutrino 
line is the dominant cosmic ray signature. Note that by gauge invariance ( II ip is 
necessarily accompanied by a factor of the hidden gauge coupling constant, while 
the operator (flOl) may not be, so that a somewhat small gauge coupling may be able 
to account for some or all of the required suppression. 

Similarly to the gravitino case, an appropriate dark matter relic abundance for 
the hidden gauge boson may be obtained through non-thermal production. For 
example, we may suppose that the inflaton decays with an appropriate branching 
fraction into the hidden sector, while also reheating the visible sector. 

4.3 A singlet fermion in an extra dimension 

Here we point out that in the context of an extra dimension, it is straightforward 
to obtain a highly suppressed coupling such as that needed for a monochromatic 
neutrino line from dark matter decay. In particular, we may consider a scenario to 
produce the operator of case 2 in table [TJ 

Suppose that there exists an Si/Z 2 orbifolded fifth dimension separating two 
branes. One of these branes, at y — 0, hosts all of the standard model fields, while 
the other, at y — £, hosts a Majorana mass term for the right handed part of a 
singlet Dirac fermion \l/ which propagates in the bulk. In addition, \I/ has a mass 
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term in the bulk and a Yukawa coupling to a lepton and the Higgs on the standard 
model brane. The zero mode of the right handed part of this bulk fermion may then 



be exponentially suppressed on the standard model brane, taking the formF 7 



*£W) = \j e 2nT_ l J M -f V ^R D \ x ) = ee m ^(x), (12) 

where M* is the fundamental scale related to the four- dimensional Planck scale by 
= M^i. Here we have written the action for the zero mode of \l/ as 

S = J d 4 xdy{M* (i¥ 0) T A d A ^ (0) + m¥ 0) ^ 0) ) 

+ [6(£ - y)M R ^H^> + S(y)x¥^L H + h.c] } . (13) 

It is then straightforward to choose Mr to be 2.4 PeV to explain the energies of the 
IceCube events, while due to exponential suppression e may be taken to be of order 
10 -29 to yield an appropriate ip dark matter lifetime even if A is of order 1. 

As in the previous examples we have discussed, inflaton decays may yield an 
appropriate dark matter relic abundance. In this case there is also another interesting 
possibility- namely, we may take the dark matter particle to carry gauged B — L 
symmetry (along with the standard model fermions, and two more right handed 
neutrino-like states for anomaly cancellation), so that the mass Mr is only produced 
after spontaneous B — L breaking. In this case, B — L interactions in the early 
universe may be used to produce the needed ip relic density. The correct abundance 
of dark matter can be attained if the reheat temperature Tr is below the scale at 
which the U(1)#-l gauge symmetry is restored, and also below the temperature at 
which ijj particles would come into thermal equilibrium through gauge interactions. 
The population of ijj particles can be produced in processes // — > ipip mediated by the 
heavy \]{1)b-l gauge boson. The resulting density to entropy ratio can be estimated 
as in Ref. 



T=T R 

3 



b 45 9* 1 

V 10 2 / Vl0 18 GeV; V5xl0 13 GeV/ ' 1 } 

where H is the Hubble parameter, g* is the number of relativistic degrees of freedom 
at the time of reheating, (av) ~ T 2 /M B _ L is the production cross section, rif ~ T 3 is 



# 7 The zero mode for the left handed part of ^ is set to zero by choosing it to be odd under the 
Z2 orbifold symmetry as usual. 
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the number density of standard model fermions in the plasma, and the first equality 
is evaluated at reheating. Numerical solution of the Boltzmann equation gives a 
value consistent with this result (13]. The dark matter mass density is then 

^0.2x( 24x 7 o * GeV ) ( TJ ^ F5 ). (15) 

5 Discussion and Future Prospects 

In this paper we have catalogued all of the operators which may lead to decays 
of PeV dark matter particles into monochromatic neutrino lines, and have given 
examples of models which may lead to appropriate decay rates to explain the two 
anomalous events recently reported by the IceCube collaboration. Here we would 
like to highlight an interesting feature of our analysis: For all of the operators that 
we have discussed in this paper one actually obtains also a lower energy continuum 
of cosmic ray neutrinos in addition to the monochromatic neutrino line. These 
are produced since in every case there are necessarily secondary decay products or 
alternate decay modes which lead to Higgses, W-bosons, Z-bosons or charged leptons 
whose decays produce a continuum of neutrino events. 

For example, in the hidden gauge boson model discussed in section H~2l the vector 
dark matter particle decays into all standard model particles carrying hypercharge. 
In particular, we will obtain decays to muons and tau leptons leading to a continuum 
neutrino signature. In the gravitino model of section 14.11 and the singlet fermion 
model in section 14.31 there are necessarily decays to W-boson + charged lepton 
which produce continuum neutrinos, in addition to those produced from the Higgs 
and Z-boson final state particles in the primary neutrino decays. For all cases we 
have considered in this paper, these final states leading to continuum neutrino signals 
have a similar branching fraction to the monochromatic neutrino events which have 
been our primary interest. We therefore have the important result that if the IceCube 
PeV events are due to dark matter decays, then there should also be a continuum of 
excess lower energy events that can also be discovered in the sub-PeV region. 

While the precise size and shape of this continuum is model-dependent, qualita- 
tively it always has a similar form. In figure [1] we show both line and continuum 
signals assuming that the partial decay width of the continuum signal is twice that 
of the line signal. This corresponds to the cases of either the gravitino model or the 
singlet fermion model discussed in the text. The combined atmospheric neutrino 
background (including those from prompt decays) [2] is also shown for comparison. 
The continuum flux was calculated using the method adopted in reference jl6], and is 
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Neutrino Energy (GeV) 

Figure 1: Line and continuum neutrino signals from PeV dark matter decays. 

based on the contribution from hadronic cascade decays of SM particles. In addition, 
we can also expect another contribution from leptonic decays, but this is not included 
in the figure for simplicity Note that in both the gravitino and singlet fermion cases 
we also have direct decays into a W-boson plus a charged lepton Z, with the flavor 
of the lepton being model dependent. Error introduced by our approximation of 
dropping leptonic decays will be negligible for the cases of I = e or I = r, while if 
I = /i, the continuum spectrum in the sub-PeV region will be somewhat enhanced. 

Let us note that there are some special cases in which the prediction of the 
continuum neutrinos may be avoided. One possibility is that the dark matter decay 
to neutrinos is taking place via an operator as in the hidden gauge boson model, 
but with this operator arising via some other means than a mixing of the new 
gauge boson with hypercharge. If the dark matter particle then couples primarily to 
electrons/electron neutrinos, and not to other flavors, then no continuum neutrinos 
result. This may seem to be a rather unlikely situation. Another possibility is that 
the dark matter particle may decay to a neutrino plus some other hidden sector 
particle. For an interesting example, one may consider a scalar dark matter particle 
which decays into two hidden sector singlet fermions ip. If ip actually mixes with 
standard model neutrinos, then this will lead to decays of <fi to ip plus a neutrino, 
without an additional continuum neutrino spectrum. We will give details of a split 
seesaw model which realizes this scenario in appendix A. 

One might wonder in addition about the possibility of other types of cosmic ray 
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signatures from decay products in our models, such as gamma rays or antiprotons. 
Unfortunately these are unlikely to be detectable in the foreseeable future. The 
reason is the following: First, backgrounds of diffuse gamma rays and cosmic ray 
antiprotons have fluxes whose energy spectra are softer than 1/E 2 , due to their 
production by cosmic-ray protons. On the other hand, gamma rays and antiprotons 
from dark matter decays have fluxes whose energy spectra are harder than 1/E 2 
(typically going as 1/E). This is because the signal spectra are essentially determined 
by the fragmentation functions of dark matter decays and these must be harder than 
1/E 2 , otherwise their integrals over energy will diverge. As a result, the ratio of the 
signal flux to the background flux becomes smaller at smaller energies. Moreover, 
both gamma rays and antiprotons are now observed at most up to 1 TeV in energy, 
making detection difficult. This situation may be clearly seen in reference [17] for 
the gamma-ray case, where it was shown that near future gamma-ray observations 
can cover dark matter lifetimes at most up to 10 27 seconds. 
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A A Split Seesaw Model 

Here we will discuss a model which is outside of the main line of argument in the 
text for two reasons: the first is that there is an additional hidden sector particle in 
the decay final state, and the second is that the decay may not be thought of as due 
to a single effective operator, since it is a result of a mixing between two low mass 
particles. As mentioned in the discussion section, the basic idea is to have a scalar 
dark matter particle (f> decaying to two light hidden sector fermions i/j through a 
(highly suppressed) (fripip interaction, and also require that ifi has some mixing with 
standard model neutrinos. We will now show that such a situation may be obtained 
in a split seesaw framework in an extra dimension jH], in which the fermion ip can 
literally be a right handed neutrino in the sense that it leads to a seesaw neutrino 
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mass in the standard model jlQl UH H2], even though ip itself will be very light P 8 

The basic setup is similar to the one used in section 14.31 We again put standard 
model fields on a y = brane in an extra dimension, with a \1> field propagating 
in the bulk as in that section, and with a zero mode wavefunction peaked on the 
brane at y — t. Again we also put a Yukawa coupling between \l/ L and H on the 
standard model brane leading to an interaction eXipLH, where we are continuing to 
use the notation of section FOl A difference here however, is that we will now put the 
Majorana mass Mr for \l/ on the standard model brane rather than the y = I brane. 
As a result, ip will obtain a highly suppressed mass of e 2 Mr. An interesting result- 
and the original motivation for the split seesaw framework- is that a seesaw mass is 
then obtained for a standard model neutrino which is interestingly independent of 
the wavefunction suppression factor e, with m v = X 2 v 2 /Mr. 

Finally, we introduce a new scalar field living on the standard model brane, 
with a Yukawa coupling to \l/ resulting in an interaction of size ge 2 (pijjip. will be 
our dark matter particle, and thus we choose its mass to be 2.4 PeV. 

After electroweak symmetry breaking, one obtains a potentially large mixing 
between ip and a neutrino v of where v is the standard model Higgs vev. While 
the primary decay mode of will be to two ip particles, as a result of the mixing, 
may also decay to ipv with a lifetime of order y l0 g ^ 8 j x 10 28 s, where we have 
taken the mixing angle to be of order 1. Obtaining an appropriate neutrino mass 
m v with A also of order 1 requires Mr to be of order 10 16 GeV as usual. Finally 



let us point out that, as was discussed in section 14. 3[ an interesting possibility for 
producing the dark matter abundance results if one assumes that carries gauged 
U(1)b-l charge, so that high temperature B — L interactions produce the relic 
particles. The estimate for the resulting relic density is analogous to that in section 



3 



# 8 We of course need more than one non-zero neutrino mass in the standard model sector, and 
thus require more than one right handed neutrino. This will not concern us here as a single ip field 
is sufficient for our present purpose. 

# 9 In the original split-seesaw papers ifi was a dark matter candidate with a keV mass which 
could explain pulsar kicks fJU 09]. However, the keV mass scale was not a definitive prediction 
of the model, which essentially creates a "democracy of scales": Majorana masses of different 
orders of magnitude can arise from the exponential suppression employed in this model. For 
example, two degenerate right-handed neutrinos with GeV masses as used in the z/MSM [50] can 
be accommodated in this scenario as well. 
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